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Abstract — The nature of Re-platinum-group element (PGE; Pt, Pd, Ir, Os, Ru) transport in the marine 
environment was investigated by means of marine sediments at and across the Cretaceous-Tertiary boundary 
(KTB) at two hemipelagic sites in Europe and two pelagic sites in the North and South Pacific. A traverse 
across the KTB in the South Pacific pelagic clay core found elevated levels of Re, Pt, Ir, Os, and Ru, each of 
which is approximately symmetrically distributed over a distance of — 1.8 m across the KTB. The Re-PGE 
abundance patterns are fractionated from chondritic relative abundances: Ru, Pt, Pd, and Re contents are 
slightly subchondritic relative to Ir, and Os is depleted by —95% relative to chondritic Ir proportions. A similar 
depletion in Os (—90%) was found in a sample of the pelagic KTB in the North Pacific, but it is enriched in 
Ru, Pt, Pd, and Re relative to Ir. The two hemipelagic KTB clays have near-chondritic abundance patterns. 
The — 1.8-m-wide Re-PGE peak in the pelagic South Pacific section cannot be reconciled with the fallout of 
a single impactor, indicating that postdepositional redistribution has occurred. The elemental profiles appear 
to fit diffusion profiles, although bioturbation could have also played a role. If diffusion had occurred over 
—65 Ma, the effective diffusivities are —10' 13 cm 2 /s, much smaller than that of soluble cations in pore waters 
(— 10 -6 cm 2 /s). The coupling of Re and the PGEs during redistribution indicates that postdepositional 
processes did not significantly fractionate their relative abundances. If redistribution was caused by diffusion, 
then the effective diffusivities are the same. Fractionation of Os from Ir during the KTB interval must 
therefore have occurred during aqueous transport in the marine environment. Distinctly subchondritic Os/Ir 
ratios throughout the Cenozoic in the South Pacific core further suggest that fractionation of Os from Ir in the 
marine environment is a general process throughout geologic time because most of the inputs of Os and Ir into 
the ocean have Os/Ir ratios >1. Mass balance calculations show that Os and Re burial fluxes in pelagic 
sediments account for only a small fraction of the riverine Os (<10%) and Re (<0.1%) inputs into the oceans. 
In contrast, burial of Ir in pelagic sediments is similar to the riverine Ir input, indicating that pelagic sediments 
are a much larger repository for Ir than for Os and Re. If all of the missing Os and Re is assumed to reside 
in anoxic sediments in oceanic margins, the calculated burial fluxes in anoxic sediments are similar to 
observed burial fluxes. However, putting all of the missing Os and Re into estuarine sediments would require 
high concentrations to balance the riverine input and would also fail to explain the depletion of Os at pelagic 
KTB sites, where at most —25% of the K-T impactor’ s Os could have passed through estuaries. If Os is 
preferentially sequestered in anoxic marine environments, it follows that the OsAr ratio of pelagic sediments 
should be sensitive to changes in the rates of anoxic sediment deposition. There is thus a clear fractionation 
of Os and Re from Ir in precipitation out of sea water in pelagic sections. Accordingly, it is inferred here that 
Re and Os are removed from sea water in anoxic marine depositional regimes. Copyright © 2003 Elsevier 
Science Ltd 


1. INTRODUCTION 

The platinum-group elements (PGE) are Pt, Pd, Ir, Os, Rh, 
and Ru. Along with Re, these elements are gaining importance 
in studies of weathering, diagenesis, and paleoclimate. For 
example, the Re-Os radiogenic decay system has been used to 
track secular changes in the relative proportions of continental, 
hydrothermal, and extraterrestrial input of Os into the oceans 
(e.g., Peucker-Ehrenbrink, 1996; Peucker-Ehrenbrink and 
Ravizza, 2000). In addition, the variable speciations of Re and 
the PGEs make them potentially useful indicators of ocean 
oxidation states (Anbar et al., 1996; Crusius et al., 1996; 
Morford and Emerson, 1999). However, interpreting Os isoto- 
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pic ratios and PGE elemental abundance data is complicated by 
an incomplete understanding of the behaviors of Re-PGE-Au 
during weathering, di agenesis, and transport in the marine 
environment. Furthermore, the full sources of Os (including 
cosmic dust infall and mantle additions from spreading centers) 
and the sinks for Os are not readily quantifiable (e.g., Sharma 
and Wasserburg, 1997; Sharma et al., 1999, 2000). 

Sediments from the Cretaceous-Tertiary boundary (KTB) 
illustrate this point. The KTB coincides with a large Ir anom- 
aly, which is now widely agreed to represent the fallout of a 
large meteorite impact. This leads to a simple prediction that 
the Ir anomaly should be manifested as a spike in the strati- 
graphic record corresponding to the time of impact and that the 
PGE abundance pattern of the spike should match that of the 
impactor, which, on the basis of the uniform relative abun- 
dances of PGEs in meteorites (Larimer and Wasson, 1988) and 
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Fig. 1. Locations of the sites analyzed in this study. Star represents K-T impact site. 


the Cr isotopic composition of the KTB sediments (Shukoly- 
ukov and Lugmair, 1998), was probably chondritic. However, 
a number of studies have shown that the Ir anomaly is often 
spread out over substantial distances in the sedimentary section 
that includes the KTB (Kastner, 1984; Kyte et al., 1985, 1993, 
1996; Robin et al., 1991; Zhou et al., 1991). In some cases, this 
distribution corresponds to several million years of deposition 
(Zhou et al., 1991), which is too long to be reconciled with a 
single impact event. There are very few KTB studies that 
consider most of the PGEs. The available data show that PGE 
abundance patterns at various KTB sites can depart signifi- 
cantly from chondritic values (e.g., Kyte et al., 1985). This 
suggests that the PGE fallout was fractionated from its initially 
chondritic proportions and redistributed before, during, or after 
deposition, or some combination of these. The goal of the 
present study is to use the distribution of Re-PGEs around the 
KTB to study the processes that fractionate and transport Re 
and the PGEs in the marine environment. The results of the 
KTB study will be combined with an investigation of Re-PGE 
content in pelagic sediments away from the KTB to draw more 
general inferences about the mechanisms of Re and PGE trans- 
port in the marine environment. 

2. SAMPLES 

2.1. Pelagic Sites 

We examined the distribution of the PGEs (Pt, Pd, Ir, Ru, 
Os) and Re at selected stratigraphic sections from KTB sites 
(Fig. 1). We analyzed a pelagic sediment core from the South 
Pacific (DSDP 596) at closely spaced intervals, including the 
KTB, which was previously identified at 20.10 m below the 
sediment-seawater interface (mbsf) on the basis of the highest 
Ir concentration peak (Zhou et al., 1991). The Ir anomaly in this 
core was found to span ~ 1 .8 m. The upper 40 m of cores at Site 
596 were recovered by hydraulic piston coring, which is con- 


sidered to have incurred minimal core disturbance. The core is 
characterized by metalliferous clays with very slow sediment 
accumulation rates (—0.3 m/Ma). All DSDP 596 sediments 
were dried in an oven at 60°C and then gently disaggregated 
and powdered by hand with a hammer covered with paper so 
that the samples always remained in a metal-free environment. 
The bulk sample was used for analyses. 

We also analyzed the KTB clay from a pelagic North Pacific 
site (DSDP 465A) dominated by calcareous oozes. The KTB 
was recognized by the shipboard scientific party (Thiede et al., 
1992) on the basis of its foraminiferal biostratigraphy 
(Boersma, 1981), and it is located in sections 3 and 4 of core 3 
in Hole 465A (—62.4 mbsf). Unfortunately, the Ir anomaly is 
smeared out over more than 0.4 m as a result of disturbance by 
rotary drilling (Kastner, 1984; Kyte et al., 1996). However, 
throughout this disturbed zone are fragments of a thin (—3 
mm), black pyrite-rich layer, and one fragment was obtained 
from 3-3 123-124. On the basis of detailed paleontological 
analyses of an intact fragment of the KTB, Kyte et al. (1996) 
showed that the black layer corresponds to the actual biostrati- 
graphic KTB. We obtained a sample, which contained the black 
layer. This boundary clay was carefully separated from the 
surrounding calcareous oozes and gently disaggregated. At the 
time of sampling, this core was nearly dry, so 465A samples 
were not oven heated. 

2.2. Hemipelagic Sites 

Two European hemipelagic KTB clays, now exposed on land 
(Caravaca, Spain, and Stevns Klint, Denmark), were examined. 
The location of the KTB and Ir peak at both of these sites has 
been previously identified (Alvarez et al., 1980; Smit and 
Hertogen, 1980). The Caravaca site is dominated by calcareous 
marls (Kyte et al., 1985). The base of the boundary clay is 
composed of a thin sanidine spherule-rich layer, in which the 
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Fig. 2. Dry weight concentrations (ppb) of Re-PGEs vs. depth 
(meters below sediment-water interface) in DSDP 596, a pelagic clay 
section in the South Pacific. Duplicate analyses are shown. Gray 
vertical line shows the location of the KTB. 


PGE anomaly is located (Smit and Hertogen, 1980). Kyte et al. 
(1985) carefully separated fragments of this boundary clay 
(SM503) from the surrounding marls and analyzed them for Re, 
Ir, Pt, Pd, and Os by radiochemical neutron activation analysis. 
In this study, we also analyzed only the sanidine spherule-rich 
boundary clay. We did not analyze the same aliquot of sample 
from Kyte et al. (1985), but we attempted to repeat the same 
sampling procedure. The separated boundary clays were com- 
bined (yielding a total of —0.5 g) and then gently disaggregated 
and homogenized by hand. The boundary clays were already 
dry, so they were not heated. 

At Stevns Klint, we analyzed a thin 1-cm layer located at the 
base of the t4 fish clay” that marks the KTB (Kastner, 1984; Kyte 
et al., 1985). The aliquot we analyzed is from the same bulk 
sample (SK10) analyzed by Kyte et al. (1985), but we also 
crushed, powdered, and homogenized our sample in a ceramic 


Spex mill (—30 g). SK10 was not heated because it was already 
dry. 


3. METHODS 

Although several PGE studies have been conducted on the above 
KTB sites, we sought to improve on these studies by analyzing Pt, Pd, 
Ir, Ru, Os, and Re on the same aliquot of each sample and by 
conducting a PGE traverse across the KTB in DSDP 596. Sample 
powders were loaded in chilled, thick-walled borosilicate glass vessels 
(Carius Tubes; Shirey and Walker, 1995) along with a mixed Re-PGE 
tracer solution ( l91 Ir, l98 Pt, ,()5 Pd, "Ru, 190 Os, l85 Re) and 3 mL of aqua 
regia (glass vessels were cleaned in boiling aqua regia). Because of the 
limited amounts of sample available, less than 200 mg were used for 
most of the dissolutions (although for samples away from the KTB, 
larger amounts were dissolved, up to 0.7 g). The glass vessel was then 
sealed and heated to 230°C for —72 h. After opening the vessel, Os was 
extracted by a three-stage chloroform solvent extraction procedure and 
then stabilized in concentrated HBr overnight (during the initial stages 
of this study. Os was not extracted and hence was lost by volatiliza- 
tion). The HBr solution was then evaporated to near dryness and taken 
up in 0.15 N HC1 for analysis of Os by inductively coupled plasma 
mass spectrometry (ICP-MS). The residual aqua regia, containing Pt, 
Pd, Ir, Ru, and Re, was then dried down and taken up in —0.5 mL of 
0.07 N HC1 and passed through 0.5 mL of cation exchange resin 
(Bio-Rad AG50X-8 100 to 200 mesh cleaned with 6.7 N HC1) followed 
by elution with 1 mL of 0.07 N HC1; the majority of interfering cations 
were adsorbed onto the resin while Re and the PGEs passed through the 
column and were collected. Approximately half the elutant was then 
used directly for ICP-MS analysis of Re, Ir, Ru, and Pt. The other half 
of the elutant from the cation exchange column was then passed 
through anion exchange chemistry (BioRad AG1X-8 100 to 200 mesh) 
to further separate Zr and Cd from Pd. A combination of 0.07 N HC1 
and 0.8 N HN0 3 washes were used to elute Zr and Cd from the column. 
Concentrated HC1 and HN0 3 were then used to elute the Re, Ir, Ru, Pt, 
and Pd, which were collected and taken up in 0.8 N HN0 3 for ICP-MS 
analysis. Analyses were performed by a magnetic-sector single-collec- 
tor ICP-MS (ThermoFinnigan Element). Procedural blanks are 3 ± 2 
pg for Os, 1 ± 2 pg of Ru, 24 ± 8 pg for Ir, 330 ± 132 pg for Pt, and 
6 ± 4 pg for Pd (2a; n = 3). The high Pt blank is likely derived from 
the glass vessel, but because of the high Pt concentrations in the 
samples, these blank levels are negligible. Over the period of this study, 
reproducibility on a 0.5 ppb mixed PGE standard of normal isotopic 
composition was better than 3 % (2a) for the isotope ratios used for 
isotope dilution calculations ( l01 Ru/"Ru, l()2 Ru/"Ru, 106 Pd/ lw5 Pd, 
l08 Pd/ l05 Pd, I87 Re/ 185 Re, l89 Os/ l90 Os, I92 Os/ 19 °Os, l93 Ir/ 191 Ir, ,94 Pt/ 
i98p t? 195 Pt/ 198 Pt, 196 Pt/ 198 Pt). Mass fractionation factors were deter- 
mined by running standards of normal isotopic composition (approxi- 
mately every 30 min) and then applying this correction factor to the 
sample measurements. Signal intensities were corrected for interfer- 
ences from Hg on 198 Pt (<1%), Cd on ,08 Pd and 106 Pd (<10%), and 
ZrO on 108 Pd and ,06 Pd (<I%). In some analyses, the Cd and/or Zr 
interference on Pd isotopes was too high (> 10%); such analyses are not 
reported. NiAr interferences on Ru were negligible, as was Pt on 192 0s. 
Spikes were calibrated against single-element Alfa-Aesar/Johnson- 
Matthey Specpure PGE standard solutions and a gravimetric Os stan- 
dard made from high purity Os salts (Lee et al., 2001; Yin et al., 2001). 
Reproducibility on different aliquots of the same sample are shown in 
Table 1. 

In a previous study of the PGE content of ultramafic rocks (Lee, 
2002), we found that the data obtained by the fire assay method agreed 
with our determinations by acid digestions in sealed glass vessels. 
However, for reasons we do not understand, we found that the fire 
assay method, when applied to DSDP 596 samples, gave poor sample- 
spike equilibration, resulting in some instances to a two- to threefold 
underestimate of the PGE concentrations compared with acid digestion 
in sealed glass vessels. 
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Fig. 3. Chondrite-normalized (Cl) Re-PGE abundance patterns at the peak (20.09 mbsf) for each KTB site. (A) 
Abundances at the peak of the Re-PGE profile at the pelagic DSDP 596 (solid lines and solid diamonds). The integrated 
concentrations minus the baseline concentrations (normalized to Cl and the measured peak at 20.09 mbsf) are denoted by 
the dotted lines and square symbols. (B) Elemental abundances of the KTB clay in the pelagic DSDP 465A core. (C) basal 
hemipelagic KTB clay from Caravaca, Spain (SM503). (D) hemipelagic KTB clay at Stevns Klint, Denmark (SK10). Values 
for chondrite normalization are Ru = 712, Pd = 560, Re -36.5, Os = 486, Ir = 481, Pt = 990 ppb (Anders and Grevesse, 
1989). Horizontal dashed lines are drawn through Ir in each panel to aid visualization. 


4. RESULTS 

4.1. Pelagic Sites 

4.1.1. DSDP 596 

Dry weight concentrations for DSDP 596 are shown in Table 
1 located at 20.09 mbsf (Fig. 2). Re, Pt, Ir, Os, and Ru are 
elevated at the KTB (Fig. 2), with the highest levels of these 
elements occurring in sample 3-4 49-50, located at 20.09 mbsf. 
Like Zhou et al. (1991), we assign the KTB to the position of 
the highest Ir concentrations. Pt, Ir, Os, and Ru levels decrease 
monotonically both up and down section from the peak, pro- 
ducing roughly symmetric concentration profiles over a width 
of —1.8 m for Pt, Ir, Os and Ru. With the exception of one 
anomalous sample (0.77 m above the peak). Re is also charac- 
terized by a roughly symmetric profile spanning —1.8 m. Most 
of our Ir measurements are in agreement to within 20% of the 
instrumental neutron activation analyses of Zhou and Kyte 
(1992). The Re, Ir, Ru, and Os contents at the KTB (20.09 
mbsf) are approximately 10, 30, 15, and 20 times higher, 
respectively, than the baseline (pre and post KTB sediments). 
Pt exhibits a slightly subdued peak, which is only a factor of 
two times higher than that of the pre-KTB sediments. There is 


a high baseline of Pt content (—25 ppb) above and below the 
KTB, and it rises to very high levels —2.5 m above the KTB. 
Pd differs from Re and the other PGEs in not displaying a 
well-defined peak at the KTB. Like Pt, the background Pd 
concentrations appear to be high (15 to 25 ppb) throughout the 
section. The rise in Pd —3 m above the KTB is coincident with 
the rise in Pt and small rises in Re, Ir, Os, and Ru. However, the 
elevated Pd contents —3.5 m below the KTB are not accom- 
panied by a rise in Re or any of the other PGEs. 

Figure 3A shows the chondrite-normalized concentrations 
(Anders and Grevesse, 1989) for the measured peak values 
(20.09 mbsf) and for the values obtained by integrating the 
concentrations over ±0.9 m around the peak after subtracting 
out the baseline concentration. For Pt and Pd, there is a signif- 
icant baseline (25 ppb for Pt and 1 5 ppb for Pd) that must be 
subtracted to obtain their contents at the peak. Ir, Pt and Pd 
have roughly chondritic proportions at 20.09 mbsf (Pt/Ir = 1 .4 
X (Pt/Ir) cr ; Pd/Ir = 0.9 X (Pd/Ir) CJ , whereas Ru, Re and Os 
appear to be significantly depleted relative to Ir (Ru/Ir = 0.5 X 
(Ru/Ir) CI ; Re/Ir = 0.6 X (Re/Ir) cl ; Os/Ir = 0.06 X (Os/Ir) CI ). 
Subchondritic Ru/Ir, Re/Ir, and Os/Ir ratios at the peak (20.09 
mbsf) are also obtained when the data are integrated over the 
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Table 1 . Dry weight concentrations. 


Ir Os Ru Pt Pd Re 

Site mbsf (ppb) (ppb) (ppb) (ppb) (ppb) (ppb) Os/Ir Pt/Ir 


Pelagic 


South Pacific - DSDP 596 
DSDP 596, 1-1 143-145 CT1 1.43 

DSDP 596, 1-2 52-54 CT1 2.02 

DSDP 596, 2-5 24-25 CT1 11 .74 

DSDP 596, 3-1 121-122 CT1 16.31 

DSDP 596, 3-2 49-50 CT1 17.09 

DSDP 596, 3-2 100-102 CT1 17.6 

DSDP 596, 3-3 122-124 CT1 19.32 

DSDP 596, 3-3 143-144 CT1 19.53 

DSDP 596, 3-3 143-144 CT2 19.53 

DSDP 596, 3-4 10-11 CT1 19.7 

DSDP 596, 3-4 29-30 CT1 19.89 

DSDP 596, 3-4 49-50 CT2 20.09 

DSDP 596, 3-4 49-50 CT3 20.09 

DSDP 596, 3-4 68-69 CT1 20.28 

DSDP 596, 3-4 68-69 CT2 20.28 

DSDP 596, 3-4 88-89 CT1 20.48 

DSDP 596, 3-4 127-128 CT1 20.87 

DSDP 596, 3-4 127-128 CT2 20.87 

DSDP 596, 3-6 100-102 CT1 23.6 

North Pacific - DSDP 465A 
DSDP 465A 3-3 123-124 CT1 62.4 


Stevns Klint, Denmark 
SK10 CT1 

SK10 (Kyte et al. 1985) 
Caravaca, Spain 
SM503 CT1 

SM503 (Kyte et al. 1985) 


0.324 

0.107 


9.45 

0.269 

0.115 


9.63 

1.48 

0.147 

1.62 

48.4 

1.61 

0.134 

1.37 

78.7 

0.427 

0.032 

1.19 

41.2 

0.769 

0.077 

0.630 

39.7 

2.57 


2.91 

26.1 

3.93 


4.01 

26.0 

4.11 

0.291 

3.64 

25.7 

6.13 


5.86 

29.1 

8.46 


7.09 

39.9 

13.6 

0.809 

10.1 

42.8 

13.8 


10.2 

38.9 

7.81 


6.15 

42.3 

7.68 

0.523 

5.85 

31.7 

5.45 


5.12 

31.0 

2.22 

0.169 


26.0 

2.20 

0.154 

1.96 

25.9 

0.489 

0.033 

0.325 

22.1 

12.7 

1.08 

27.2 

117 


Hemipelagic 


50.9 

32.2 

62.8 

76.4 

47.4 

35 


78 

42.0 

31.6 

53.7 

85.9 

56.9 

46 


131 



0.156 

0.331 

29.2 


0.232 

0.426 

35.8 

4.16 

0.112 

0.099 

32.7 


0.181 

0.083 

48.9 

22.7 

0.057 

0.075 

96.6 


0.135 

0.100 

51.6 


0.325 


10.1 

17.2 

0.247 


6.62 

11.7 

0.206 

0.071 

6.25 

22.3 

0.313 


4.75 

19.1 

0.480 


4.71 

19.7 

0.704 

0.059 

3.15 

19.4 

0.623 


2.82 

16.2 

0.354 


5.42 

17.4 

0.357 

0.068 

4.12 

19.7 

0.286 


5.69 

16.8 

0.181 

0.076 

11.7 

15.8 

0.176 

0.070 

11.8 

24.0 

0.112 

0.067 

45.2 

90.1 

9.41 

0.085 

9.18 


44.7 

8.30 

0.633 

1.5 

46 

9.4 

0.738 

1.6 

72.4 

1.80 

0.752 

2.0 

10 

2.7 

0.808 

2.3 


a Concentrations are based on dry weight; mbsf = meters below sediment-water interface; duplicates were not run in all cases; SM503 represents 
a thin (~1 mm) sanidine-spherule-rich boundary clay. 


width of the Re-PGE anomaly (—1.8 m). However, because of 
the high baseline levels for Pt and Pd, the integrated Pt/Ir (0.7 
X (Pt/Ir) CI ) and Pd/Ir ratios (0.4 X (Pd/Ir) CI ) are subchondritic, 
indicating that the apparent chondritic Pt/Ir and Pd/Ir ratios at 
the peak may be fortuitous. 

Up-section from the KTB, the Re-PGE relative abundances 
are also fractionated from chondritic values. Pt/Ir (14 to 47 X 
(Pt/Ir) CI ) and Pd/Ir ratios (2 to 42 X (Pd/Ir) CI ) are distinctly 
elevated compared with chondrite. Os/Ir ratios remain subchon- 
dritic throughout the upper part of the section rising from 0.08 
X (Os/Ir) CI at — 20 m below the sediment-seawater interface to 
0.43 X (Os/Ir) CI at ~2 m below the sediment-seawater inter- 
face (Fig. 4). The rise in PGE abundance, particularly in Pt 
abundance, at 16.31 mbsf coincides exactly with a rise in cobalt 
content, as determined by Zhou and Kyte (1992). 

4.12. DSDP 465 A 

Table 1 shows the dry weight Re-PGE concentrations deter- 
mined for the thin reduced layer marking the KTB at the 
pelagic DSDP 465 A. Our concentration determinations are 
consistent with those of Kyte et al. (1996). Chondrite-normal- 
ized values are plotted in Figure 3B. Relative to chondritic 
proportions, Pt, Pd, and Re are highly enriched relative to Ir 
(Pt/Ir = 3.2 X (Pt/Ir) CI ; Pd/Ir = 4.6 X (Pd/Ir) CI ; Re/Ir = 8.4 X 
(Re/Ir) CI ). The Ru/Ir ratio is only slightly superchondritic (1.4 
X (Ru/Ir) CI ). These large enrichments in Pt, Pd, and Re, and 


moderate enrichment in Ru, contrast with the lack of large 
enrichments seen in the integrated KTB composition at DSDP 
596. As we did not measure samples above and below the KTB, 
we are unable to estimate the magnitude of the “baseline” 
contribution of Pt, Pd, and Re to the KTB at DSDP 465A. The 
one feature in common between 465 A and 596 is their dis- 
tinctly subchondritic Os/Ir ratios; the Os/Ir ratio of the reduced 
layer at 465A is 0.09 X (Os/Ir) CI . Considering the high Os and 
Ir contents, we consider the Os/Ir ratio of the 465A boundary 
clay to be robust. 

4.2. Hemipelagic Sites 

4.2.1. Caravaca 

The Re-PGE concentrations (Table 1 and Fig. 3C) obtained 
for sample SM503 are ~33% lower than that of Kyte et al. 
(1985) on similar samples, but the relative abundances are in 
very good agreement. The discrepancy is most likely due to 
sample heterogeneity. Caravaca differs from the pelagic DSDP 
sites in having roughly chondritic PGE relative abundances 
(Pt/Ir = 0.9 X (Pt/Ir ) CI , Pd/Ir = 1.1 X (Pd/Ir) CI , Ru/Ir = 0.8 X 
(Ru/Ir) CI ; Os/Ir = 0.8 X (Os/Ir) CI ) and a subchondritic Re/Ir 
ratio (0.5 X (Re/Ir) CI ). Although we did not measure above and 
below the boundary clay, Kyte et al. (1985) reported declining 
PGE abundances in the 3 cm immediately above SM503 and 
low concentrations in background samples above and below the 
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Fig. 4. Os/Ir ratios from this study (circles) vs. mbsf in DSDP 596 core. Squares represent the combined Ir and Os 
measurements of Zhou and Kyte (1992) and Peucker-Ehrenbrink (1996), respectively (thick horizontal line = chondritic 
Os/Ir ratio). Upper horizontal axis provides approximate age (Ma) of sediment (Zhou and Kyte, 1992). Curves represent the 
Os/Ir ratio in pelagic sediments as a function of a change in the area fraction (F) of anoxic sediment deposition. From bottom 
to top, each curve represents a linear decrease from F = 0.2% at the present sediment-seawater interface to F = 4.0, 1.5, 
0.8, 0.4, and 0.2%, respectively, at 25 mbsf. These calculations are based on Eqn. 6a, with 33 pg/m 2 yr used as the value 
for the pelagic Os burial flux (Table 4), 1.5 X 10 5 pg/m 2 yr for the anoxic Os burial flux, and the riverine Os/Ir ratio given 
in Table 2. 



KTB. They found that Os/Ir ratios actually decline somewhat in 
PGE-rich sediments above SM503, but only to ratios —50% 
lower than in chondrites. They did not observe the extremely 
low Os/Ir ratios seen here in pelagic sites. 


4.2.2. Stevns Klint 

Re-PGE abundances for the boundary clay from the hemipe- 
lagic Stevns Klint section are shown in Table 1 and chondrite- 
normalized abundances are shown in Figure 4d. PGE abun- 
dances agree to within 10% and Re abundances to within 15% 
of Kyte et al. (1985). Considering that the samples were pro- 
cessed separately, this is remarkable agreement. Except for Re, 
the PGE abundances have near-chondritic proportions (Pt/Ir = 
1.5 X (Pt/Ir) CI , Pd/Ir = 0.9 X (Pd/Ir) cl , Ru/Ir = 1.2 X 
(Ru/Ir) cl ; Os/Ir = 0.7 X (Os/Ir) CI ), and therefore Stevns Klint 
resembles Caravaca. Re is somewhat elevated with respect to Ir 
(1.9 X (Re/Ir) CI ). We did not measure above and below the 
boundary clay, but given the high peak concentration, it is 
likely that any baseline levels of Re-PGEs are comparably 
small. 


5. DISCUSSION 

We first discuss the origin of the broad (—1.8 m) Re, Pt, Ir, 
Os, and Ru anomaly at DSDP 596 (Fig. 2). This depth interval 
corresponds to an apparent deposition time of —6 Ma (based on 
a constant sedimentation rate), which is much too long to be 
reconciled with the fallout of a single impact event and there- 
fore requires postdepositional Re-PGE redistribution in the 
stratigraphic column. We will also discuss the origin of the high 
degree of Re-PGE fractionation in the two pelagic sites com- 
pared with the two hemipelagic sites, in particular, focusing on 
the origin of the extremely low Os/Ir ratios in the two pelagic 
sites. 


5.1. Causes of Postdepositional Re-PGE Mobilization at 
DSDP 596 


One way to obtain a symmetric profile for Re, Pt, Ir, Os, and 
Ru is by diffusional redistribution of an initial spike in PGE, 
mathematically represented by 


C = 


— t=Z= ] e -iA 7(4Df) 
p J4 nDt ) 


( 1 ) 


where C is the concentration (after subtraction of the back- 
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from 2 X JO 2 to 2 X 10 4 yr, for Os ranges from 1 X 10 5 to 6 
X 10 6 yr, and for Re ranges from 2 X 10 8 to 1 X 10 9 yr. The 
large discrepancies between residence times calculated from 
the riverine input (t 2A ) and the pelagic burial output ( t 2B ) 
reflects the fact that Os and Re are primarily being sequestered 
in other sinks besides the deep sea. The high Os and Re r 2B s, 
however, appear to be geologically unreasonable, suggesting 
that oceanic margins, rather than estuaries, represent the other 
sink for Os and Re. The similarity between r 2B and r 2A for Ir 
(—2 X 10 3 yr; Table 2) is consistent with pelagic sediments 
representing the main sink for Ir input into the ocean. 

Finally, there is also the possibility that Os and Re can be 
removed from seawater when it returns to the estuary (Martin 
et al., 2001). In the extreme case, when only seawater Os is 
removed and no riverine Os is removed, the oceanic residence 
times would remain unchanged from scenario A. 

6. CONCLUSIONS 

The nature of Re-PGE transport in the marine environment 
was investigated by using marine sediments at and across the 
KTB at two hemipelagic and two pelagic sites. A traverse 
across the KTB in a Pacific pelagic clay core found elevated 
levels of Re, Pt, Ir, Os, and Ru, each of which is approximately 
symmetrically distributed over a distance of — 1.8 m across the 
KTB. The 1.8-m interval containing the Re-PGE anomaly 
cannot be reconciled with the fallout of a single impactor with 
chondritic abundances, indicating that postdepositional redis- 
tribution has occurred. If the elemental profiles are the result of 
diffusion that occurred over —65 Ma, then the effective diffu- 
sivities are —10“ 13 cm 2 /s, much smaller than that of soluble 
cations in pore waters (— 1CK 6 cm 2 /s). Regardless of the cause 
of redistribution, the fact that Re and the PGEs are spread out 
over the same lengthscale indicates that postdepositional pro- 
cesses did not significantly fractionate their relative abun- 
dances. If redistribution was caused by diffusion, then effective 
diffusivities of Re and the PGEs must be equal. The Re-PGE 
abundance patterns are fractionated from chondritic relative 
abundances, both at the KTB and integrated over the whole 
section. Os is depleted by —95% relative to chondritic Ir at the 
pelagic KTB in the South Pacific. A similar depletion in Os was 
found in a sample from another pelagic KTB section in the 
North Pacific. However, the two hemipelagic KTB clays have 
near-chondritic abundance patterns. Because the fractionations 
seen in the pelagic sections predate deposition, fractionation of 
Os from Ir during the KTB interval must have occurred during 
aqueous transport in the marine environment. Extremely sub- 
chon dri tic Os/Ir ratios throughout the Cenozoic in the pelagic 
South Pacific section further suggest that fractionation of Os 
from Ir in the marine environment is a general process through- 
out geologic time because most of the inputs of Os and Ir into 
the ocean have Os/Ir ratios ^ 1 . 

The generality of this process is confirmed by a global mass 
balance calculation, which shows that pelagic Os and Re burial 
fluxes account for only a small fraction of the Os (<10%) and 
Re (<0.1%) riverine inputs into the ocean while Ir burial in 
pelagic sediments appears to be a major sink for Ir in the 
oceans. If all of the missing Os and Re is assumed to reside in 
anoxic sediments in oceanic margins instead of pelagic sedi- 
ments, then the respective burial fluxes calculated by mass 


balance for these sediments is comparable to observed values 
of Os and Re in organic-rich sediments in oceanic margins. In 
contrast, putting all of the missing Os and Re into estuarine 
sediments would require high concentrations to balance the 
riverine input and would also fail to explain the depletion of Os 
at pelagic KTB sites, where at most —25% (the area of land 
surface) of the K-T impactor’ s Os could have passed through 
estuaries. There is thus a clear fractionation of Os and Re from 
Ir in precipitation out of sea water in pelagic sections. Accord- 
ingly, it is inferred here that Re and Os are removed from sea 
water in anoxic marine depositional regimes. If Os is prefer- 
entially sequestered in anoxic marine environments, then the 
Os/Ir ratio of pelagic sediments should be sensitive to changes 
in the rates and areas of anoxic sediment deposition. 
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(6a) 

Eqn. 6a depends on the Os/Ir ratio of the riverine input rather 
than the absolute Os and Ir concentrations in rivers. Taking 0 Os 
and 0 ,r to be —0.2 (Sharma et al., 1999) and the removal rate 
of Ir by anoxic sediments is negligible (e.g., p 3 a 2 3^j r / 
P4CK24C4O <<1), Eqn. 6a simplifies to 

C? /C\ (1 - F) 

c'l ~ V C'J ( Pl a n C?/p 4 a 24 C?) F + ( 1 - F) (6b) 

Thus, when the area fraction of anoxic sediments is negligible 
(F = 0), the Os/Ir ratio of pelagic sediments is equal to the 
riverine Os/Ir ratio. When F approaches 1, the Os/Ir ratio of 
pelagic sediments approaches zero. Because of the strong af- 
finity of Os for anoxic sediments, only a very small amount of 
anoxic sediment deposition is needed to decrease the Os/Ir ratio 
of pelagic sediments. It follows that a shift in Os/Ir ratio of 0.06 
to 0.4 in the upper part of the DSDP 596 section may simply 
reflect a decrease in deposition of anoxic sediments. We have 
illustrated in Figure 3 how a decrease with age in the area of 
anoxic sediment deposition, F, from —1.5% at —80 Ma to 0.2% 
today can reproduce the pattern of Os/Ir ratios in the DSDP 596 
section, for the case in which p 3 a 23 Cf y / p 4 a 24 C 4 s and the 
Os/Ir ratio of the rivers remain constant and j3°7/3 ,r — 1. 
Because Re also has a high affinity for anoxic marine sedi- 
ments, variations in rates of anoxic/suboxic sediment deposi- 
tion should lead to covariations in Re/Ir and Os/Ir ratios. Re/Ir 
ratios in DSDP 596 indeed show an increase with decreasing 
depth below the sediment-seawater interface. It is possible that 
variations in anoxic sediment deposition may also control the 
shifts in other PGEs seen in DSDP 596. 

If the Os and Ir in the upper part of the DSDP 596 core are 
not directly derived from seawater but are instead derived from 
a continental detrital component with Os/Ir > 1 , the increase in 
Os/Ir ratios from <0.1 to —0.4 could be due to an increase in 
detrital proportion rather than reflecting changes in seawater 
composition. As the concentration of Ir in the upper —2 m of 
the DSDP 596 section is —300 pg/g (Table 1), a detrital 
component corresponding to — 1 00 pg/g Ir and — 1 00 pg/g Os 
is needed to increase the Os/Ir ratio of the sediment from 0. 1 to 
0.4. These values are much higher than the Os and Ir concen- 
trations of upper continental crust (—22 pg/g Ir and —31 pg/g 
Os; Peucker-Ehrenbrink and Jahn, 2001). Thus, although the 
sedimentation rates in the DSDP 596 core appear to have 
increased with younging age (Zhou and Kyte, 1992), the above 
calculation shows that the Os/Ir ratios in these sediments pre- 
dominantly reflect seawater and not the detrital component. 

5.6. Further Implications 

Understanding the transport properties of Re and the PGEs in 
the marine and estuarine environments has important implica- 
tions for calculating their residence times in the ocean, that is 
the average amount of time the element of interest resides in the 
ocean. For the case in which there is no removal of an element 
in estuaries (scenario A; Fig. 7 A), the inverse of the element’s 
oceanic residence time, t 2A , can be expressed as 


where M 2 is the mass of the ocean, C 2 is the concentration of 
the element of interest in the ocean, and all other symbols are 
previously defined. Eqn. 7 shows that the inverse of the oceanic 
residence time, for the case in which removal is completely 
within the oceans, is equal to the summation of the rate con- 
stants for deposition in the oceanic sinks or to the rate constant 
of the riverine input. In most cases, the residence time is 
calculated from the latter because the total of all the oceanic 
sinks is not well known. Assuming 0 is small (0 — 0.2; Sharma 
et al., 1999), estimates of the riverine Re-PGE input and the 
Re-PGE concentration in seawater (Table 2) can be used to 
calculate in this manner. This yields 2 X 10 3 , 4 X 10 4 , and 
7 X 10 3 yr for the oceanic residence time of Ir, Os, and Re, 
respectively. These values are in accord with previous studies 
in which estuarine removal was also considered to be small or 
zero (Anbar et al, 1992, 1996; Colodner et al., 1993; Sharma 
and Wasserburg, 1997; Levasseur et al., 1999; Sharma et al., 
1999). If the dominant removal process of Ir, Os, and Re is 
within the oceans instead of the estuaries, then the burial flux in 
oceanic margins Q 24 C 4 can be calculated from the known 
pelagic burial and riverine input fluxes. 

On the other hand, if removal occurs in the estuaries and in 
only one oceanic sink (scenario B; Fig. 7B), the residence time 
is fundamentally changed. The inverse of the oceanic residence 
time t 2b is expressed as 


1 

7 2 B 
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"/.3+f » 


CoiCo 
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QiaC* 
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where f 13 is the fractional removal of riverine flow in the 
estuary and all other symbols are previously defined. Because 
Q 0 ,C 0 in scenario B is equivalent to Q 02 C 0 in scenario A, Eqn. 
8 can be reexpressed to yield the scenario B oceanic residence 
time, r 2/? , in terms of the scenario A residence time, r 2A : 


_ (1 + fi)T 2 A 
Tw ~ (1 + J3) ~fn 


(9a) 


Because 0 is small, Eqn. 9a simplifies to the approximation 


T 2 «~t^ (9b) 

Thus, if estuarine removal occurs, the residence time calculated 
assuming transport processes described by scenario A (Eqn. 7) 
must be increased by a factor of 1/(1 - f 13 ). For example, if 
>90% and >99.9% of the riverine Os and Re, respectively, are 
sequestered in estuaries, scenario B residence times for Os and 
Re will be orders of magnitude higher, i.e., 5 X 10 5 and 7 X 
10 6 yr (scenario B) instead of 5 X 10 4 and 7 X 10 3 yr (scenario 
A), respectively. As the fraction of riverine component re- 
moved in estuaries, f l3 , is not known very well, scenario B 
oceanic residence times (r 2B ) can be calculated from pelagic 
burial fluxes (right-hand side of Eqn. 8) because these sedi- 
ments represent the only oceanic sink in scenario B. Scenario B 
residence times (r 2B ) for Ir, Os, and Re calculated in this 
manner are shown in Table 3 along with the pelagic burial flux 
estimates used in the calculations. We find that r 2B for Ir ranges 
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at least 1 .4 X 1 0 7 pg/m 2 yr. These oceanic margin burial fluxes 
for Os and Re are roughly in agreement with direct estimates of 
the Os (6 to 12 X 10 4 pg/m 2 yr) and Re (1.6 X 10 7 pg/m 2 yr) 
burial fluxes in organic-rich sediments found in oceanic mar- 
gins (Ravizza and Turekian, 1992; Colodner et al., 1993), 
consistent with anoxic marine sediments being a significant 
repository for Os and Re (Ravizza and Turekian, 1992). Su- 
perchondritic Os/Ir ratios measured on black shales (Os/Ir — 4 
to 21 for the Ordovician Utica black shales in Quebec; Peucker- 
Ehrenbrink and Hannigan, 2000), are also consistent with Os 
being preferentially sequestered over Ir in anoxic environ- 
ments. 

If scenario A is the dominant mechanism for redistributing 
Os and Re in the ocean (i.e., over 90% of oceanic Os and 
essentially all of the Re is transported from the open ocean to 
oceanic margins, which make up only a small fraction of the 
entire ocean floor), the physical plausibility of this scenario 
must be assessed. A minimum constraint can be placed on the 
mass flow of water in the oceans. Assume that a wall is placed 
in the ocean for scenario A, such that the deep open ocean is 
designated by reservoir 2 and oceanic margins, where Os is 
sequestered (reservoir 3), is designated by reservoir 1 (Fig. 7C). 
The mass flow of water between reservoirs 1 and 2, Q , 2 and 
Q 21 , are equal to conserve mass (assuming no differences in 
evaporation rate from reservoirs 1 and 2). The mass balance of 
Os for the open ocean (C 2 ) is 0 = ( 1 -f j3)Q 02 C 0 - Q 24 C 4 + 
Q , 2 (C , - C 2 ). The pelagic burial flux of Os, Q 24 C 4 , is very 
small, so this equation simplifies to Q 12 = (1 + /3) (Q 02 C 0 )/ 
(C 2 -Cj). A minimum constraint on the mass flow in the ocean, 
Q 12 , is given by Q l2 — Q 02 Q) (1 + /3) /C 2 . In the case of the 
modern oceans and rivers, the concentrations of Os in rivers 
and the open ocean are nearly equal (C 0 = C 2 ), and therefore, 
Q 12 — Q 02 (1 + j3). Thus, to move nearly all of the Os from the 
deep ocean to another region in the ocean for deposition, the 
water flow in the ocean must be equal to or greater than the 
global riverine water input (Q 02 ), which is 3.6 X 10 16 kg/yr. If 
it is assumed that the concentration of Os in the oceans varies 
by no more than —10% (e.g., (C 2 - Cj)/C 0 — 0.1)), then Q )2 
= 3.6 X 10 ' 7 kg/yr for the case in which fi = 0 (it is unlikely 
that ft is greater than 1 ; Sharma et al., 1999). Taking the mixing 
time of water in the oceans to be — 1500 yr (Broecker and Li, 
1970), the average mass flow of water in the oceans is therefore 
estimated by dividing the mass of water in the oceans by the 
mixing time. This gives a mass flow of water in the oceans of 
9 X 10 17 kg/yr. These calculations show that redistribution of 
Os completely within the oceans by ocean currents is physi- 
cally plausible. The intraocean current flow rates indicated here 
should be sufficient to redistribute other elements in the open 
ocean to localized regions where they could be sequestered. 

We now assess Os transport in scenario B where all of the 
missing Os is sequestered in estuaries instead of oceanic mar- 
gins. If the area covered by estuaries is —3.8 X 10 n m 2 , or 
0.1% of the area of the oceans (Woodwell et al., 1973), we can 
use Eqn. 3, taking /3 = 0, to calculate the minimum Os 
concentration of estuarine sediments given some estimate of 
the sediment mass accumulation rate. Martin et al. (2000) 
measured the Os concentration of estuarine sediments at the 
mouth of the Fly and Sepik rivers, which drain the island of 
Papua New Guinea. If we generalize the sediment mass accu- 
mulation rates reported by Martin et al. (2000), 1 .5 to 15 X 10 3 


g/m 2 yr, we find that estuarine sediments should have at least 50 
to 500 pg/g of Os. These minimum values are slightly higher 
than the estuarine sediment concentrations (86 ± 43 pg/g; 1 
SD) reported by Martin et al. (2000), suggesting that estuarine 
sediments may not be the dominant sink for Os. The implau- 
sibility of scenario B can also be revealed by examining Eqn. 
3. In the extreme case in which riverine Os is completely 
removed in estuaries, f , 3 = 1 , the ensuing requirement that ^24 
= j3 (Eqn. 2) would indicate that the isotopic composition of Os 
in seawater is governed solely by cosmic infall or hydrothermal 
fluxing. This is certainly not allowed as the Os isotopic com- 
position of seawater is much more radiogenic than meteoritic 
and hydrothermal signatures (Sharma and Wasserburg, 1997; 
Levasseur et al., 1998; Sharma et al., 1999). These conclusions 
are further corroborated by the fact that Os/Ir ratios are ex- 
tremely subchondritic at the KTB in pelagic sediments. As 
oceans make up —75% of the Earth’s surface, at most 25% of 
the impactor’s Re-PGE fallout could have passed through riv- 
ers and estuaries if the fallout was homogenously distributed 
over the Earth’s surface. 

We consider that Os transport is best described by scenario 
A, in which most of the Os and virtually all of the Re is 
sequestered in oceanic margins rather than in oxic pelagic 
environments. Ir appears to be strongly sequestered in pelagic 
sediments instead of oceanic margin sediments. This implies 
that the Os content of pelagic sediments underestimates the 
extraterrestrial influx, whereas the Ir content of pelagic sedi- 
ments overestimates the extraterrestrial influx (Peucker-Ehren- 
brink, 1996). 

5.5. Significance of Variations in Os/Ir in Pelagic Clays 

One significant feature of the DSDP 596 core is that the Os/Ir 
ratio in the core sediments rises from 0.06 to 0.1 near the KTB 
to —0.4 at —2 m below the sediment-water interface in (Fig. 
4). The change in Os/Ir ratio of pelagic sediments must be due 
to changes in the Os/Ir ratio of seawater. We explore the simple 
case in which the Os/lr ratio of the total oceanic input is 
constant and the changing Os/Ir ratio of the ocean is due to 
changes in the deposition rates of Os and Ir in the ocean as a 
result of changes in the areal distribution of anoxic and oxic 
sediments. Assuming steady state, the mass balance of Os in 
the ocean (Fig. 7 A) can be expressed as 

Qo 2 Co s 

(1 + p° x ) — = p 4 a 24 CfV — F) + p 3 a 23 C?F, (4) 

where A is the area of the seafloor, F is the fraction of the total 
area occupied by anoxic sediments [F = A 3 /(A 3 4- A 4 )], and all 
other terms are defined as before. Clearly, the concentration of 
Os in the sediments C c 3 )s and C 4 v will change for a fixed 
riverine input of Os (G 02 C^ V ) if the sediment deposition rates 
change. It follows that changes in the area of anoxic sedimen- 
tation would cause substantial changes in the Os concentration 
in pelagic sediments. If Eqn. 4 is rewritten for Ir, 

0 - )C !r 

(1 + pf)-~ - = p 4 a 24 <t( 1 — F) + Pi a 23 C' 3 r F, (5) 

we can combine Eqn. 4 and 5 to yield a general expression for 
the Os/Ir ratio of pelagic sediments. 
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Fig. 7, Box models describing two different scenarios in which there 
is a missing sink. (A) Scenario A. Both sinks are located in the oceans, 
where Qo 2> Q24> and Q 23 are the mass flow of river water (reservoir 1 ) 
into the ocean (reservoir 2), pelagic sediment accumulation (reservoir 
4), and sediment accumulation in hemipelagic environments (reservoir 
3). C 0 , C 2 , C 3 , and C 4 represent the concentrations of an element in 
river water, seawater, hemipelagic sediment, and pelagic sediment, 
respectively. (B) Scenario B. One of the sinks occurs between the river 
and the ocean, i.e. in estuaries (reservoir 1), thereby reducing the actual 
amount of an element entering the open ocean (reservoir 2). Variables 
as in scenario A, except Q 01 , Q 12 , and Q 13 represent the river flow into 
the estuarine water reservoir, mass accumulation rate of estuarine 
sediments, and water flow from the estuarine reservoir to the ocean, (c) 
Scenario A modified by placing a wall in the ocean, generating two 
boxes, 1 and 2- Q|2 and Q 2 , represent the mass flow of water between 
boxes 1 and 2. 


does not include the fraction of extraterrestrial influx to conti- 
nents, which is assumed to eventually enter the ocean via 
rivers. The mass balance of Os in the oceanic reservoir (reser- 
voir 2) is then given by 0 = (1 + jS)Q 02 C 0 - Q 23 C 3 - Q 24 C 4 . 
Dividing this equation by Q 02 C 0 and letting f 23 = Q 2 3C 3 /Qo2Q> 


and f 24 = Q 2 4C 4 /Q 02 C 0 , the following form of the equation can 
be obtained, where the input and output terms are expressed in 
fractional form relative to the riverine input term 

f 2 3 + f24= l+P (2) 

Eqn. 2 simply states that the fractional amount of riverine Os 
buried in oceanic margins (f 23 ) and in the deep sea (f 24 ) is equal 
to one plus the fractional amount of nonriverine input relative 
to the riverine input /3. 

In scenario B (Fig. 7B), Os is removed when rivers (reservoir 
0) pass through estuaries en route to the ocean (reservoir 1), 
thus diminishing the actual amount of Os released to the open 
ocean (reservoir 2). Reservoir 3 represents the estuarine sedi- 
ment and reservoir 4 represents the pelagic sediment. No other 
deep-sea sink is included in this scenario. As in scenario A, we 
define the fraction of all nonriverine inputs directly into oceans 
relative to the riverine input as j3 = (total mass flow of 
nonriverine inputs directly to the oceans)/(Q ol C 0 ). The Os mass 
balance for the oceanic reservoir (2) is then given by 0 = 
PQo\ C 0 + Q| 2 C] - Q 24 C 4 and that for the estuarine reservoir 
(7) is given by 0 = Q 0 iC 0 - Q 13 C 3 - Q 12 C,, where Q 13 and 
Q 24 are the estuarine and pelagic sediment mass accumulation 
rates. We assume that nonriverine input into estuaries with 
concurrent removal is negligible. The mass balance equation 
for reservoir 1 can be rearranged to give (1 — f 13 ) = (Q 12 C,/ 
Q 0 ,C 0 ) where f 13 is the mass fraction of riverine Os lost from 
the river and sequestered in estuarine sediments (reservoir 3), 
e.g., fn = (Q, 3 C 3 )/(Q 0 iC 0 ). The mass balance equation for the 
oceanic reservoir can then be expressed as 0 = Q 01 C 0 (1 - f 13 ) 
+ £Q0,C 0 — Q24^ 4 . Dividing by the unmodified riverine input 
Qo.C 0 and expressing the fractional burial of riverine Os in 
pelagic sediments as f 24 = Q24CVQ01C0 we obtain the frac- 
tional form of the oceanic reservoir mass balance: 

f, 3 + f 24 = 1+0 (0<fj 3 < 1 ) (3) 

Although Eqn. 3 appears identical in form to Eqn. 2, there is 
a subtle but important difference. Because f 13 represents the 
fractional amount of riverine Os removed in estuaries, f 13 
cannot be greater than 1 because the nonriverine inputs into the 
ocean, by construct of this model, never interact with estuaries 
(Fig. 7B). In the extreme case of f 13 = 1, there is no riverine 
input of Os into the ocean, and thus, Eqn. 3 requires f 24 0» 

that is the pelagic burial of Os equals the total nonriverine Os 
inputs. 

A minimum estimate of the Os removal rates in oceanic 
margins in scenario A can be determined from mass balance 
(Eqn. 2) by using the pelagic burial rates and riverine input 
rates given in Table 4. This calculation shows that pelagic 
sediments account for <10% of the total riverine Os input, 
<0.1% of the total Re input, and all of the Ir input to within 
error. If scenario A describes Os transport, then at least 90% of 
the Os and essentially all of the Re must be buried in oceanic 
margin environments. If anoxic marine sediments represent the 
other sink for Os and Re in the oceans, the present areal 
distribution F of anoxic sediments in the oceans (—0.3%) can 
be used to estimate the minimum burial flux of Os and Re in 
this missing sink (Table 4). The Os burial flux in anoxic 
sediments required to balance the pelagic sediment burial flux 
is found to be at least 3 X 10 5 pg/m 2 yr, whereas that for Re is 
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Fig. 6. Comparison of Re, Os, and Ir input (Table 2) and pelagic 
output fluxes for different lithologies (Table 3). The minimum input 
flux is represented by the riverine input (arrows denote lower bound). 
Each data point assumes that the entire ocean floor is covered by a 
particular lithology. The average burial of Os in deep-sea sediments, 
assuming a 1:1 ratio of calcareous oozes and pelagic clays, is —115 
pg/m 2 yr and is denoted by the open diamond. The two points for 
DSDP 596 refer to burial flux estimates for portions of the stratigraphic 
section having different sedimentation rates. Bold diagonal line repre- 
sents one-to-one correspondence, in which the pelagic burial flux 
balances the riverine input (dashed lines represent deviations from the 
one-to-one line at 1 0-fold increments). 


Ehrenbrink and Ravizza, 2000). Exactly where the “missing” 
Re and Os lies is of immediate concern because it bears directly 
on understanding how Re and Os are transported in the hydro- 
sphere. In particular, the calculated residence times of Re 
and/or Os in seawater depend critically on the number of Re 
and Os sinks and their respective removal rates. If estuarine 
removal is the dominant removal process of riverine Re and Os, 
the calculated residence times will be significantly different 
than in the case in which the removal occurs predominantly in 
a single uniform deep ocean sink. One means of identifying the 
extra sinks can be had by studying estuaries and oceanic margin 
environments themselves. However, in the case of Os, the 
relative importance of such sinks is not known very well. 
Martin et al. (2000, 2001) and Levasseur et al. (2000) show that 


dissolved Os is removed from waters in estuarine environ- 
ments. In detail, however, there are subtle, but important, 
differences between these studies. Martin et al. (2001) suggest 
that a significant component of oceanic Os may be sequestered 
in tropical estuaries, whereas Levasseur et al. (2000) argue that 
a large amount of dissolved riverine Os is sequestered in 
estuarine sediments. Thus, even though studying oceanic mar- 
gins and estuaries clearly yields valuable information, the 
highly variable depositional and chemical conditions in oceanic 
margins and estuaries makes it difficult to representatively 
sample Os deposition in these environments. In contrast, the 
spatial variation of Re-PGE concentrations in the open oceans 
is more constant. In this section, pelagic sediments and simple 
mass balance principles are therefore explored to construct a 
more generalized model for Re, Os, and Ir transport in the 
hydrosphere. 

We examine two end-member steady-state scenarios describ- 
ing the transport of Re-PGEs. The formulations are focused on 
Os, but can be applied to other elements. In scenario A, the 
extra Os sink is taken to occur in oceanic margins (e.g., 
hemipelagic environments, suboxic and anoxic marine basins, 
etc.). In scenario B, the extra Os sink is taken to occur in 
estuarine environments. Both of these scenarios are represented 
as box models, as shown in Figure 6. 

In scenario A (Fig. 7A), Os from rivers (reservoir 0) flows 
into the ocean (reservoir 2) and is sequestered in pelagic 
sediments (reservoir 4) and oceanic margin sediments (reser- 
voir 3). We also include the net contribution of other inputs of 
Os into the ocean (cosmic infall, hydrothermal influx, Aeolian 
infall, and groundwater flow), which are poorly quantified. The 
ocean is assumed to be well-mixed. The mass flow (kg/yr) of 
river water into oceans is denoted by Q 02 , and the mass accu- 
mulation rates of pelagic and oceanic margin sediments in the 
oceans are given, respectively by Q 24 and Q 23 , where for 
sediments is equal to a^-Ay (a- is the sedimentation rate from 
i to j , pj is the dry weight density of the sediment, and Aj is the 
area covered by the jth reservoir). Let C 0 , C 3 , and C 4 be the 
concentrations of the element of interest in river water, oceanic 
margin sediments, and pelagic sediments, respectively. We 
define the fraction of all nonriverine inputs into the ocean 
relative to the riverine input as f3 = (total mass flow of 
nonriverine inputs directly into the ocean)/(Q 02 C 0 ). Defined in 
this manner, nonriverine inputs represent only direct inputs into 
the ocean. For example, in the case of extraterrestrial influx, /3 


Table 4. Calculated burial fluxes in oceanic margins based on average petagic and riverine fluxes. 


Element 

Lithology 

Average burial 
flux in pelagic 
sediments 
(pg/m 2 yr) 

Riverine flux 
(pg/nryr) 

Fraction of total 
burial required in 
ocean margin 
sinks 

Burial flux in anoxic 
sediments 
(F = 0.3%) 
(pg/m 2 yr) 

Ir 

Pelagic clay 

103 

175 

0.41 

2.4 X 10 4 

Ir 

Calcaraeous ooze 

147 

175 

0.16 

9.3 X 10- 

Ir 

1:1 pelagic — calcareous 

125 

175 

0.29 

1.7 X 10 4 

Os 

Pelagic clay 

33 

898 

0.96 

2.9 X 10 5 

Os 

Calcaraeous ooze 

198 

898 

0.78 

2.3 X 10" 

Os 

1:1 pelagic — calcareous 

115 

898 

0.87 

2.6 X I0 5 

Re 

Pelagic clay 

23 

42681 

0.9995 

1.4 x 10 7 

Re 

Calcaraeous ooze 

87 

42681 

0.9980 

1.4 X 10 7 

Re 

1:1 pelagic — calcareous 

55 

42681 

0.9987 

1.4 X 10 7 




664 


C.-T. A. Lee, G. J. Wasserburg, and F. T. Kyte 


Table 3. Pelagic burial fluxes. 11 


C 4 a 24P^4 T 2B 


cone, in a 24 pelagic apparent 
Dry sed sed. rate burial flux residence 


Element and site 

Reference 


Lithology 

Location 

(pg/g) 

(m/Ma) 

pg/m 2 yr 

time (yr) 

Ir 

DSDP 596 (>5 Ma) 

This study 


Pelagic clay 

S Pacific 

400 

0.4 

61 

6 

X 

10 3 

DSDP 596 (<5 Ma) 

This study 


Pelagic clay 

S Pacific 

300 

0.4 

46 

9 

X 

10 3 

GPC 3 (15-60 Ma) 

(Kyte et al., 1993) 


Pelagic clay 

N Central Pacific 

1800 

0.2 

137 

3 

X 

10 3 

GPC 3 (<1.5 Ma) 

(Kyte et al., 1993) 


Pelagic clay 

N Central Pacific 

200 

2.2 

167 

2 

X 

10 3 

DSDP 577 (near KTB) 

(Michel et al., 1985) 


Calcareous ooze 

NW Pacific 

20 

3 

23 

2 

X 

I0 4 

DSDP 292 (30-40 Ma) 

(Keller et al., 1987) 


Calcareous ooze W Equatorial Pacific 

20 

45 

342 

1 

X 

10 3 

DSDP 462 (30-40 Ma) 

(Keller et al., 1987) 


Clay 

Equatorial Pacific 

40 

135 

2052 

2 

X 

I0 2 

DSDP 612 (30^10 Ma) 
Os 

DSDP 596 (>5 Ma) 

(Keller et al., 1987) 


Calcareous ooze 

W Atlantic 

10 

20 

76 

5 

X 

10 3 

This study 


Pelagic clay 

S Pacific 

40 

0.4 

6 

6 

X 

10 6 

DSDP 596 (<5 Ma) 

This study 


Pelagic clay 

S Pacific 

100 

0.4 

15 

3 

X 

I0 6 

GPC3 (15-60 Ma) 

(Pegram et al., 1992) 


Pelagic clay 

N Central Pacific 

120 

0.2 

9 

6 

X 

10 6 

GPC3 (<1.5 Ma) 

(Pegram et al., 1992) 


Pelagic clay 

N Central Pacific 

70 

2.2 

59 

3 

X 

10 6 

16-162-17-4, 139-140 (48.7 
Ma) 

20-196-2-5, 112-114 (80 Ma) 

(Peucker-Ehrenbrink et al.. 

1995) 

Chalk 

Equatorial Pacific 

38 

3 

43 

9 

X 

10 5 

(Peucker-Ehrenbrink et al.. 

1995) 

Pelagic clay 

NW Pacific 

182 

1.25 

86 

4 

X 

10 5 

39-2-6, 130-131 (53 Ma) 

(Peucker-Ehrenbrink et al.. 

1995) 

Pelagic clay 

NE Pacific 

416 

0.34 

54 

7 

X 

10 5 

1 29-80 1A-4R-2-45-47 (57.1 
Ma) 

9-77B-53-1, 70-72 (37.9 Ma) 

(Peucker-Ehrenbrink et al.. 

1995) 

Pelagic clay 

NW Pacific 

9 

1 

3 

1 

X 

10 7 

(Peucker-Ehrenbrink et al., 

1995) 

Pelagic clay 

E Equatorial Pacific 

64 

10 

243 

2 

X 

10 5 

34-319-12-3, 134-136 (15.8) 

(Peucker-Ehrenbrink et al., 

1995) 

Calcareous ooze 

E Equatorial Pacific 

61 

16.5 

382 

1 

X 

10 5 

92-597-6-2, 106-107 (26.3 Ma) 

(Peucker-Ehrenbrink et al., 

1995) 

Calcareous ooze 

E Equatorial Pacific 

34 

7.7 

99 

4 

X 

10 5 

8-74-12-3, 18-20 (42 Ma) 
Re 

DSDP 596 (>5 Ma) 

(Peucker-Ehrenbrink et al.. 

1995) 

Calcareous ooze 

Central Pacific 

75 

2.3 

66 

6 

X 

\0 5 

This study 


Pelagic clay 

S Pacific 

150 

0.4 

23 

1 

X 

10 9 

DSDP 596 (<5 Ma) 

This study 


Pelagic clay 

S Pacific 

150 

0.4 

23 

1 

X 

10 9 

8-74-12-3, 18-20 (42 Ma) 

(Peucker-Ehrenbrink et al.. 

1995) 

Calcareous ooze 

Central Pacific 

38 

2.3 

33 

9 

X 

10 8 

92-597-6-2, 106-107 (26.3 Ma) 

(Peucker-Ehrenbrink et al.. 

1995 

Calcareous ooze 

E Equatorial Pacific 

48 

7.7 

140 

2 

X 

10 8 


a Notation follows that given for scenario B transport model (reservoir 0 = river; reservoir 2 = ocean; reservoir 4 = pelagic sediment); C 4 = dry 
weight concentration in sediment (pg/g); a 24 = sedimentation rate (m/Ma); pelagic burial flux is Q 24 C 4 or a 24 pC 4 , where p is the dry weight density 
(assumed to be - 0.38 g/cm 3 ); r 2B = residence time for transport model of scenario B = (C 2 M 2 )/(a 24 pC 4 A 2 ), where A 2 is the area of the ocean (3.61 
X 10 14 m 2 ), M 2 is the mass of the ocean (1.4 X 10 21 kg), and C 2 is the concentration in the ocean (Table 2); sediment accumulation rates and ages 
are based on biostratigraphic data and can be found from the Initial Reports of the Deep Sea Drilling Project (volumes 8, 9, 16, 20, 34, 39, 86, 91, 
92, and 129) or from Kyte et al. (1993) and Zhou and Kyte (1992) for GPC 3 and DSDP 596, respectively. 


oozes and if we generalize the limited data set, then it can be 
seen that burial fluxes of Os and Re in calcareous oozes are 
higher, whereas that for Ir does not change to within error. The 
higher burial fluxes for Os and Re in calcareous oozes is due 
mainly to the high sedimentation rates (1 to 2 orders of mag- 
nitude higher than pelagic clay accumulation rates, Table 3). 
Presently, calcareous oozes make up <50% of the Pacific 
Ocean, —60% of the Atlantic Ocean, and —50% of the Indian 
Ocean (Davies and Gorsline, 1976). Assuming a 1:1 ratio of 
calcareous oozes to pelagic clays, the average deep-sea burial 
flux of Os would be 1 15 pg/m 2 yr, which is shown as a diamond 
in Figure 6 and tabulated in Table 4 (average pelagic clay burial 
flux for Os is 34 and that for calcareous oozes is 198). This 
value accounts for only — 1 3% of the riverine input of Os (900 
pg/m 2 yr). Given that our estimate of the Os input is a minimum 
bound for the total Os input into the oceans, Os burial in 
deep-sea sediments in the present-day oceans clearly does not 
balance the Os input. Moreover, the present day areal distribu- 
tion of pelagic calcareous sediments may be the highest it has 
ever been since the Cretaceous (Opdyke and Wilkinson, 1988). 
Opdyke and Wilkinson (1988) showed that the total rate of 
pelagic limestone accumulation is approximately eight times 
higher at the present as compared with the end of the Creta- 


ceous. Calcareous oozes therefore cannot solve the problem of 
missing Os in the deep-sea. 

In summary, we conclude that deep-sea sediments represent 
a much larger sink for Ir than for Os and Re. Mass balance 
considerations assuming the riverine Re-PGE inputs represent a 
minimum estimate of the total inputs into the ocean indicate 
that at least 95% and 99.9% of the total oceanic input of Os and 
Re, respectively, must be sequestered in sinks other than deep- 
sea sediments. These extra sinks could be located in oceanic 
margins or estuaries. The amount of Ir sequestered in these 
extra sinks must be much smaller than Os and Re, but the exact 
amount is not known because the total oceanic input of Ir is 
unconstrained. In the next section, we use simple models to 
speculate where the “missing” Os and Re resides. 

5.4. Implications for the Transport of Re, Ir, and Os 

The last section demonstrated that other sinks for Re and Os 
must exist besides pelagic sediments. Numerous studies have 
shown that Re and Os can be sequestered in oceanic margins 
(e.g., anoxic and suboxic sediments in hemipelagic environ- 
ments or marginal basins) and in estuaries (Colodner et al., 
1993; Levasseur et al., 2000; Martin, 2000, 2001; Peucker- 
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Table 2. Riverine fluxes. 3 


Element 

C„ 

Cone, in 
river (pg/kg) 

Qo2^(/A 2 

Riverine flux 
(pg/m 2 yr) 

C 2 

Cone, in 
ocean (pg/kg) 

t 2 a 

residence 
time (yr) 

Ir 

1.75 

175 

0.1 

2 X 10’ 

Os 

9 

898 

10 

4 x 10 4 

Re 

428 

42681 

8000 

7 X 10’ 


a Notation follows that used for scenario A box model (reservoir 0 = 
river; reservoir 2 - ocean); C 0 = concentration in river water; Q 02 = 
global river flow (3.6 X 10 16 kg/yr); A 2 = area of ocean (3.61 X 10 14 
m 2 ); t m = residence time for scenario A = (C 2 M 2 )/(Q 02 C 0 ), where M 2 
is the mass of the ocean (1.4 X 10 21 kg); oceanic and river concentra- 
tions of Os, Ir, and Re are from Anbar et al. (1992), Colodner et al. 
(1993), Anbar et al. (1996), Sharma and Wasserburg (1997), Levasseur 
et al. (1999), and Peucker-Ehrenbrink and Ravizza (2000). 


fractionation of the K-T impactor’s Os and Ir was probably 
widespread in the marine environment. The majority of Os 
from the K-T impactor therefore had to be soluble, allowing it 
to be separated from Ir (which could be either in a soluble or 
insoluble form) and deposited elsewhere besides the deep sea. 
This conclusion is consistent with that made by Peucker-Ehren- 
brink (1996), which is based on fewer data for DSDP 596. 

It is likely that this fractionation of Os from Ir during the K-T 
interval is also a general process throughout geologic time. The 
Os/Ir ratio throughout the rest of the DSDP 596 section is far 
less than the chondritic ratio of 1 , ranging from 0.06 at the KTB 
to 0.4 just below the sediment-water interface (Fig. 3). Sub- 
chondritic Os/Ir ratios have also been observed in other pelagic 
sediments (Koide et al., 1991; Ravizza and Pyle, 1997; Ravizza 
et al., 2001). These subchondritic ratios in pelagic sediments 
contrast with the chondritic to slightly superchondritic Os/lr 
ratios in the main inputs of Os and Ir into the oceans (extra- 
terrestrial influx M (Brownlee et al., 1984; Peucker-Ehren- 
brink, 1996), continental crust ~1 (Peucker-Ehrenbrink and 
Jahn, 2001), rivers ~3 to 16 (Anbar et al., 1996; Levasseur et 
al., 1999; Peucker-Ehrenbrink and Ravizza, 2000), there is at 
present no direct estimate of the Os/Ir ratio of the hydrothermal 
input, but the Os/Ir ratio of the mantle is ~1). Aqueous frac- 
tionation of Os from Ir is therefore spatially and temporally 
widespread (Ravizza et al., 2001). Furthermore, it can be qual- 
itatively concluded that deep-sea sediments cannot account for 
all of the Os input into the oceans. In the next section, the sinks 
and sources of Os, Ir, and Re will be quantified by means of 
simple transport models. 

5.3. Oceanic Mass Balance of Os, Ir, and Re 

In this section, we seek to further assess the generality of 
Os/Ir fractionation in the marine environment by quantifying 
the magnitudes of the deep-sea Os and Ir sinks, and comparing 
these values to the global input of Os, Ir, and Re into the 
oceans. The global input of Os, Ir, and Re into the oceans can 
come from rivers, extraterrestrial infall (in vapor or particle 
form), hydrothermal fluxing, groundwater flow or Aeolian in- 
fall. Only the riverine input is reasonably well-known, and 
therefore, we use the riverine input as a minimum estimate of 
the Os, Ir, and Re inputs into the ocean. The global riverine flux 
of Os, Ir, and Re are tabulated in Table 2, where the riverine 


flux is equal to Q riv C riv /A oc (Q riv is the global riverine mass 
flow [kg/yr], C riv is the concentration [pg/kg] of an element in 
river water, and A nc is the area of the oceans [m 2 ]). These 
riverine fluxes are derived from riverine concentrations prior to 
any estuarine removal. In Table 3, we list the bulk pelagic 
sediment burial fluxes calculated from estimated sedimentation 
rates and bulk concentrations for various localities and pelagic 
lithologies in the Pacific Ocean by using the available data. The 
average age of the pelagic sediments used for these calculations 
is ~40 Ma. The burial flux was expressed as C' sed pa, where 
C' sed is the dry weight concentration (pg/g) of element i in the 
sediment, p is the sediment dry weight density (g/cm 3 ), and a 
is the sedimentation rate (mm/yr). Sedimentation rates were 
taken from the Initial Reports of the Deep Sea Drilling Project 
(volumes 8, 9, 16, 20, 34, 39, 61, 86, 91, 95, 92, and 129) or 
from Kyte et al. (1993) and Zhou and Kyte (1992). The sedi- 
mentation rates are based on biostratigraphic markers. In Figure 
6, pelagic burial outputs vs. riverine inputs are compared on a 
logarithmic scale. Because the global riverine inputs of Os and 
Re are based on several studies encompassing most of the 
world’s largest rivers (Colodner et al., 1993; Sharma and Was- 
serburg, 1997; Levasseur et al., 1999; Peucker-Ehrenbrink and 
Ravizza, 2000). the uncertainties on the estimates of the global 
Os and Re riverine inputs are considered to be insignificant in 
a logarithmic plot. The global riverine input of Ir, however, is 
less constrained because the data set is limited to two rivers 
entering the Baltic Sea and the variation in the Ir content of 
these two rivers is significant (Anbar et al., 1996). The range in 
riverine Ir input is denoted in Figure 6. For the output fluxes, 
each of the values given in Table 3 are plotted for various 
pelagic lithologies and localities for Os, Ir, and Re. In con- 
structing Figure 6, each point reflects the deep-sea burial flux 
assuming the entire ocean is covered by a particular lithology. 
The bold diagonal line represents the condition in which the 
deep-sea burial fluxes equal the riverine input. Points that plot 
to the left of this line indicate that the pelagic burial flux for a 
particular lithology is less than the riverine flux. Adding other 
inputs into the ocean will shift the data upward (as shown by 
the arrows in Fig. 6), increasing the discrepancy between total 
oceanic inputs and burial in pelagic sediments. 

It can be seen from Figure 6 that the burial flux of Os and Re 
associated with pelagic clays is 1 to 2 and 2 to 3 orders of 
magnitude lower than the riverine influx of Os and Re, respec- 
tively, whereas the pelagic clay burial flux of Ir is equal to the 
riverine influx of Ir to within uncertainties. A mass balance 
calculation based on the average of the pelagic clay burial 
fluxes and the riverine fluxes for Os and Re (Table 4) shows 
that pelagic clays account for no more than 4% and 0.05%, 
respectively, of the total inputs of Os and Re into the ocean. For 
Ir (Fig. 6), the pelagic burial fluxes are equal to the riverine flux 
to within the uncertainty of the riverine Ir input (Anbar et al., 
1996). If most of the deep-ocean is underlain by pelagic clays, 
then the data shown in Figure 6 indicate that pelagic clays serve 
as a much larger sink for Ir than for Re and Os. This generality 
is in accordance with the observations made in the previous 
section that there is a clear deficit in Os relative to Ir in pelagic 
clays at and away from the KTB in two pelagic sections. 

We now consider calcareous oozes as a deep-sea sink for Os, 
Re, and Ir. Calcareous ooze data points are shown as triangles 
in Figure 6. If the entire ocean were covered with calcareous 



ln(Ru - 0.33) ln(lr-0.55) 


662 


C.-T. A. Lee, G. J. Wasserburg, and F. T. Kyte 





0 0.2 0.4 0.6 0.8 

(AX ) 2 

Fig. 5. Natural logarithm of the baseline-corrected concentration ln(C - C B ) vs. A* 2 , where Ax is the absolute distance 
(m) from the peak. The background concentration is estimated from the data points away from the KTB Re-PGE anomaly. 
Error bars incorporate the two standard deviation of duplicate runs and an arbitrary 10% error on our choice of the 
background correction. The diffusional redistribution model of an instantaneous flux of Re-PGEs predicts a linear 
relationship in ln(C - C B ) - Ax 2 , with a slope equivalent to - 1/(4Z>0, where D is the effective diffusivity (cm 2 /s). 
Effective diffusivities for Re, Pt, Ir, Os, and Ru are shown for a diffusional timescale of 65 Ma. For Re and Pt, linear 
regression lines were made by ignoring one and two outliers (gray diamonds), respectively. Also shown are the empirical 
fits of the data to an exponential decay function, C = C 0 e -I ^ ,/L and plotted in ln(C - C B ) - A* 2 space (curved dashed line). 
Fits to the exponential function include all data points. Exponential fits are as follows: Ir = 12.98e^ 248A * (r 2 = 0.98), Os 
= OJ79t~ 220 ^ (r 2 = 0.99), Ru = 9.86e _1 96il,r (i r 2 = 0.92), Re = 0.536e~ 249 ^ (r 2 = 0.96), and Pt = 21.4e“ 217 ^ ( r 2 
= 0.89). 
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ground), C T is the total amount of the PGE per unit area 
(ng/m 2 ), p is the dry weight density, D is the effective diffu- 
sivity (m 2 /s), t is time (s), and Ax is distance (m) away from the 
peak (Crank, 1975). Such a diffusion profile yields a linear 
relationship when InC is plotted against Ax 2 , with the slope of 
the line being inversely proportional to the diffusivity (slope = 
— MADt). Figure 5 shows that the concentration profiles yield 
strong linear correlations in InC — Ax 2 space for Ir, Ru, Os, and 
Re, and weak correlations for Pt. We did not attempt to plot Pd 
in InC - Ax 2 space because of the large uncertainty in calcu- 
lating the background-corrected Pd concentration. The strong 
linear correlations for Ir, Ru, Os, Re, and to a lesser extent Pt 
are suggestive of diffusional redistribution. Assuming diffusion 
has operated since the time of deposition (—65 Ma), the effec- 
tive diffusivities are calculated to be between 2.3 to 5.3 X 
10“ 13 cm 2 /s for Re, Pt, Ir, Ru, and Os. These effective diffu- 
sivities are significantly lower than those of soluble cations in 
pore waters (10~ 6 cm 2 /s; McDuff, 1981), implying Re-PGE 
diffusion is limited by chemical reaction (limited solubility, 
adsorption, etc.) between the solid phases and interstitial pore 
waters. In making these calculations, we eliminated one Re and 
two Pt data outliers (Fig. 5). However, including these points 
does not change the above range of calculated diffusivities, i.e., 
the diffusivities remain identical to within a factor of about two. 
If the redistribution of Re-PGE at the KTB was indeed con- 
trolled by diffusion, the similarity in effective diffusivities 
indicates that Re-PGE abundances cannot be significantly frac- 
tionated at the peak because elemental fractionation at the peak 
is inversely proportional to the square root of the ratio of the 
elemental diffusivities. 

Although the diffusion model appears to describe the data, 
Zhou et al. (1991) report the presence of shocked quartz in one 
spot -0.5 m below the Ir peak, which clearly requires mechan- 
ical redistribution, such as bioturbation. However, bioturbation 
of an initial PGE spike predicts an exponentially decreasing 
PGE profile (Berger and Heath, 1968), which is inconsistent 
with the symmetric Re-PGE profile. Bioturbation should lead to 
exponentially decreasing Ir concentrations with increasing dis- 
tance up-section, e.g., C = C () e _r/L , where z is the distance (or 
the sedimentation rate X time) from the base of the initial 
homogenization zone of width L (located a distance L below 
the initial deposition layer), and C 0 is the concentration at z = 
0. We note that both the down-section and up-section legs of 
the Re, Pt, Ir, Os, and Ru profiles can each be fit to an 
exponentially decaying function (dashed lines in Fig. 5), but 
such profiles have no physical relevance to the simple biotur- 
bation model described above, which predicts one-sided expo- 
nential decay, and hence an asymmetric profile. In addition, the 
calculated bioturbation lengthscale is —0.5 m, which is much 
larger than the 0.02- to 0.1 -m lengthscale observed in most 
sediments (Thomson et al., 2000). If bioturbation is responsible 
for postdepositional redistribution of Re-PGEs, a more com- 
plex bioturbation model is clearly needed to explain the data. 
Besides the one sample lying —0.5 m below the Re-PGE peak, 
Zhou et al. (1991) did not examine other points away from the 
immediate vicinity of the peak for shocked quartz. It is there- 
fore not possible to determine if this particular sample contain- 
ing shocked quartz indicates widespread homogenization by 
bioturbation or if this sample fortuitously encountered an iso- 


lated worm hole, which extends well below the homogenization 
zone (Lowemark and Werner, 2001). 

Regardless of what caused the redistribution of Re-PGEs, the 
fact that these elements were redistributed over the same 
lengthscale indicates that they remained remarkably coupled 
during redistribution. Thus, the fractionated Re-PGE composi- 
tions were not imposed by postdepositional redistribution pro- 
cesses. Instead, fractionation of Re-PGE at the KTB appears to 
have occurred in the aqueous marine environment. 

5.2. Fractionations of Os/Ir at Pelagic KTB Sites 

The most striking result of this study is that the two pelagic 
KTB sediments have highly fractionated Re-PGE patterns, 
contrasting with the two hemipelagic KTB sediments, which 
have near-chondritic PGE patterns. The largest and most con- 
sistent fractionation appears to be a depletion in Os relative to 
Re and the other PGEs, as exemplified by extremely subchon- 
dritic Os/Ir ratios (0.06 to 0.08) for the KTB at DSDP 596 and 
465A. This result agrees with the combined Os measurements 
by Peucker-Ehrenbrink et al. (1995) and the Ir measurements 
by Kyte et al. (1996). Relative to the chondritic value of Os/Ir 
= 1, the measured Os/Ir ratios correspond to 90 to 95% 
depletion in Os. Similar magnitude depletions in Os relative to 
Ir have been seen at the KTB in a pelagic Caribbean Sea section 
(ODP Site 1001; Pearson et al., 1999) and a central Pacific 
pelagic site (LL44-GPC3; Kyte et al, 1993; Pegram and 
Turekian, 1999). Collectively, this suggests that subchondritic 
Os/Ir appears to be a widespread feature of deep-sea KTB 
sediments. 

One possibility for nonchondritic Os/Ir is that the K-T im- 
pactor had a nonchondritic Os/Ir ratio. This seems unlikely 
because the Cr isotopic composition of KTB sediments is 
similar to that of carbonaceous chondrites (Shukolyukov and 
Lugmair, 1998), and most meteorites (chondrites and metals) 
have an Os/Ir ratio of —1 (Larimer and Wasson, 1988). In 
addition, this hypothesis would require an alternative explana- 
tion for the chondritic PGE abundance pattern of the two 
hemipelagic sites analyzed here. If we assume the K-T impac- 
tor had chondritic Re-PGE abundance patterns, then the frac- 
tionations must have occurred in the marine environment. In 
the case of DSDP 596, postdepositional redistribution of Os in 
the sedimentary column cannot explain the deficit in Os relative 
to Ir at the KTB because there is not enough Os anywhere in the 
core to balance the deficit. To illustrate, the Os/Ir ratio calcu- 
lated by integrating over the elevated Os and Ir profiles (—1.8 
m) associated with the KTB remains strongly subchondritic. In 
addition, even if all of the Os in the entire stratigraphic column 
is integrated and combined with only the total amount of Ir 
present in the — 1.8-m-wide interval containing the KTB (total 
Os in sedimentary column calculated assuming a total sediment 
thickness of 70.5 m and a constant Os concentration of 0.1 
ppb), the Os/Ir ratio is raised from 0.06 to only 0.33, which is 
well below the chondritic value. 

The above calculations require the “missing” Os from the 
KTB in DSDP 596 to be deposited elsewhere besides this 
deep-sea sediment core and without Ir. Combined with the fact 
that the only four pelagic KTB sections analyzed so far for both 
Os and Ir are also characterized by distinctly subchondritic 
Os/Ir ratios (DSDP 465 and 596, LL44-GPC3, ODP-1001), 



